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In his renormalization paper of the Kondo effect Wilson replaced 
the full band of s-electrons by a small number of "Wilson states". He 
started from a rather artificial symmetric band with constant density 
of states and constant interaction with the impurity. It is shown in 
the present paper that with a minor modification the Wilson states 
are optimally suited to treat the interaction of an impurity with an 
arbitrary s-band. Each Wilson state represents electrons of a whole 
energy range Ae. It carries the interaction of all these electrons with 
the impurity. All the other electron states in this energy range have 
zero interaction with the impurity and are neglected in the calcula- 
tion. The resulting error is minor. As an example the singlet-triplet 
excitation energy of a Kondo impurity is numerically calculated for a 
tight-binding band with a strongly energy dependent density of states. 
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1 Introduction 



The description of metal electrons as plane waves ipk (r) = e^^^/ ^/V has been 
extremely successful during the past. However, if one wants to manipulate 
the whole electron basis one runs into problems. For example, a rotation 
of the whole basis in Hilbert space requires a number of Euler angles which 
is equal to N [N — 1) /2 where N is the number of basis states. Such a 
rotation is rather difficult for = 10^^ in a macroscopic metal sample. 
Wilson [1] encountered a similar problem when he treated the Kondo effect 
where the impurity mixes the plane waves in complex ways. For his numerical 
evaluation Wilson strongly reduced the number of states to a few states and 
only included these "Wilson" -states in his treatment. (Wilson called these 
states "Kondo" -states). 

The purpose of this paper is to illuminate some of the physics behind the 
use of the Wilson states and to introduce a generalization. We will consider 
a host with a d-impurity. Wilson used a band of electrons with constant 
density of states and a constant interaction with the impurity. This may 
appear to be rather artificial and far removed from a real system. This is far 
from true. With a minor modification the Wilson states can be applied to 
bands with arbitrary energy density of states and impurity interaction. This 
generalization does not increase the complexity of the calculation. 

To demonstrate the method two examples will be treated quantitatively 
in this paper. Our group developed recently a new and very compact approxi- 
mate solution for the Friedel- Anderson and the Kondo impurity problems [2] , 
[?] , [3] . This method yields remarkably accurate results for the ground-state 
energy and the singlet-triplet excitation energy. The ground state requires 
the construction of one localized s-electron state (per spin) aj, = ^ckq while 
the remaining electron states are rearranged into states at so that the new 
basis is orthonormal. Since this requires an orthogonal transformation 

of the whole electron basis Wilson's state reduction is essential in the numer- 
ical procedure. In section IV the ground-state energy and the singlet-triplet 
excitation energy for a Kondo impurity is numerically calculated for (i) a 
free electron system with p (e) oc ^/6 and (ii) a tight-binding metal. For both 
electron systems the density of states is not constant. Particularly in the 
tight binding band p{e) show a strong variation (see Fig.2). 
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2 The Wilson States 



Let us first start (like Wilson) with the somewhat artificial energy band hav- 
ing the half-band width D, i.e., the range [—D : D] with a constant density 
of states. It is half filled and the Fermi level hes at zero energy. We divide all 
energies by the half-width D of the band (see Fig.l). Then the band extends 
from —1 to +1. The full band can accommodate one electron per atomic 
volume (and spin). The density of states (per spin) is po = 1/2 pei' atomic 
volume va- 

Wilson sub-divided the positive and negative parts of the energy band into 
N cells. This is done on an exponential scale. The boundaries of the cells 
are for u < N/2 at = -1/2", i.e. -1, -1/2, .. - l/2\ .... and for u > N/2 
at = +1/2^-", i.e. ...,1/22,1/2, 1. (Wilson used originally il/A'^ with 
the parameter A but generally chose the value A = 2). With these discrete 
points one obtains a sequence of intervals or cells. The interval u (for 
iy<N/2) is given by ^.-i = -1/2^^-1 <s < -\/2" = or [^.-i : U (The 
equivalent intervals for positive values are treated completely analogously). 

The new Wilson states ipi^ lie in the center of the energy interval {^u-i, ^u) 
and have an energy of £^ = {$,^ + ^^-i) /2, i.e. -|, -|, ••, 
(for u < N/2). One essential property of the Wilson basis is that it has an 
arbitrarily fine energy spacing at the Fermi energy. In his normalization 
treatment of the Kondo effect Wilson entangled these states numerical with 
the Kondo impurity. 

The first cell €i extends from (-1 : -1/2) . This cell contains Zi(w lO^^) 
electron states which we denote as (/?^ (r). They all interact with the impurity. 
Now these states are transformed with the matrix A where A has the simple 
form Axfj, — exp {i2nXfj,/Zi) /-sfZ^. With A\^j, one obtains the new states 



where < A < Zi and the sum is taken over all states ^p^^ (r) in d. 



XA (r) = ^E^M (r) e^''^'^/^'^ 
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Fig. la: Left side: Wilson's rectangular symmetric energy band. 

It is divided into energy cells whose width decrease 

exponentially towards the Fermi energy. Each energy cell v 

is represented by a state whose energy lies in the center 

of the cell. It represents the full interaction of all states in the cell 

with the impurity. 

Fig. lb: A realistic band with energy dependent density of states. 
The division into energy cell is completely analogous. The 
construction of the Wilson states is modified to accommodate 
the energy dependence of density of state and s-d-interaction. 

If all states <y9^ (r) have the same s-d matrix element Vsd with the impurity 
then only the state for A = has a finite interaction with the impurity. For 
all other states with A 7^ the s-d interaction is zero. These states are 
ignored in the calculation. Then the cell Ci is represented by a single state 
T/^i (r) where 

^i(r) = E^^(r)/v^ (1) 
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The energy (expectation value) of this state follows from the energies of 
the states </7^ 

This state has the s-d-matrix element Vgd (1) = \fZ\Vsd- 

The same procedure is then apphed to all the other cells C,^. In other 
words each energy cell is represented by just one state which carries the 
full interaction with the impurity. The neglected states have no interaction 
with the impurity. 

For the discussed case (constant density of state and constant interaction) 
the s-d-matrix elements Vgd iy) of the Wilson states are proportional to \fZl 
and given by 

v;.(-) = K°.V^ 

where A^,^ = ^^—^v-i is the width of the energy cell and |K?;|^ = |Kd f 

I |2 

The size of the energy cells is about equal to the distance from the Fermi 
level. From perturbation theory one realizes that each Wilson state has about 
the same coupling to the impurity. (The ratio of the square of the matrix 
element divided by the energy denominator (z/)|^ I yields about the 
same contribution as long as !\ey can be replaced by itself.) 

Wilson's philosophy was somewhat different when he introduced the "Kondo"- 
states. In his renormalization approach he considered the wave functions of 
these states as spherical shells, like in an onion, surrounding the impurity. 
The states ■0,^,'?/'Ar+i-j/ belonged to the [y — l)-th shell. During the renormal- 
ization procedure he first coupled the states ipi^ipp^ directly to the impurity. 
Then he progressed step wise outwards by coupling the impurity complex 
(consisting of the impurity plus {u — 1) shells) to the shell v. In our prag- 
matic application we consider the Wilson states as representatives of a finite 
energy range which contain the complete s-d-interaction of all the (neglected) 
states in the energy range with the impurity. The square of the resulting 
Vsd {i^) is proportional to the total number of states in that energy interval. 

The neglected states which do not interact with the impurity cause a 
small error since they interact through the kinetic energy with the states 
if) I,. But Wilson pointed out that this error is very small, a result that our 
calculations confirm. 
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3 Generalization of Wilson States 



The Wilson concept can now be applied to an s-band with arbitrary density 
of states and state-dependent s-d-interaction. Numerically it is a trivial ex- 
tension of the Wilson rectangular band. Let us consider a band that extends 
from —Di to +D2 with the Fermi energy at zero energy. The density of 
states depends strongly on the energy and the matrix elements for the s-d- 
interaction Vsd are not only a function of the energy but depend on the exact 
quantum numbers of the state. To simply the discussion we normalize the 
energy with D — max{Di, D2) so that the band hes within [—1 : -|-1] (see 
fig. lb). Then we subdivide the band as before obtaining the same energy 
cells CCj,. Again we use the cell <ti to demonstrate the new procedure. This 
cell contains Zi electron states ip^ (r). 

3.1 Friedel- Anderson impurity 

Let us first consider the case where the electrons interact with the impurity 
through the s-d-interaction. Now this interaction Vsd,ij. with the impurity 
depends on the state (/?^. In this case we define the Wilson state ijji (r) as 



(The sum over /i is taken over all states (p^ within the energy cell Ci.) The 
energy ei of this state ipi follows from the energies C/j, of the states (p^ 



This state has the full s-d-interaction of the cell fti with impurity. All 
other {Zi — 1) states which can be constructed from the states 99 ^within cell 
Ci (being orthonormal to each other and to ipi (r)) have zero interaction 
with the impurity. This can be easily seen with the following argument. The 



state ipi has the coefficients a^' = Vsd,iJ./yYl I'^sd./il basis of {v?^t}- If 

X = ^P^^n is a state orthogonal to ijji (built from the same basis {^/j,}) 
then the scalar product has to be zero, i.e. 





£1 = 



/ y \'^sd,n\ 



(3) 





1 
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The s-d-matrix element between the state x ^^nd the impurity is equal to the 
sum of its components, i.e. 

From the orthogonality with respect to the state ipi follows the disappearance 
of the interaction. Therefore all other states built from the {y^^j-basis have 
no interaction with the impurity. 

So now we have a new Wilson basis {il'u} with the energies and the s-d- 
interaction Vsd {v)- The energies don't lie anymore in the center of the cells 
and the s-d-matrix elements are no longer proportional to the square root of 
the cell width. (There may be even a cell or two with no Wilson state in it. 
Then one can either reduce the number of Wilson states or one gives these 
states zero s-d- interaction.) For a numerical evaluation this generalization 
of the Wilson states to arbitrary density of states and interaction causes no 
additional complications. 

3.2 Kondo impurity 

For a Kondo impurity the situation is slightly more difficult because the inter- 
action matrix element (yj^ (r) | J (r)| ^p^i (r)) is taken between two s-electrons. 
We restrict ourselves here to the case where the exchange potential is a 5- 
function J (r) = J5(r). In this case the matrix element takes the value 
Jip*^ (0) ip^i (0) and can be separated. Therefore the role of Vsd,^l in equ.(I2]) is 
replaced by ip^ (0). Then the Wilson state ipi for the Kondo impurity is 

^1 (r) = E^M (0) V, (r) / ^JyMM? (4) 

These Wilson states have a maximal amplitude at the origin. All other 
states in £i vanish at the origin and therefore they don't interact with the 
5-exchange potential. 

The matrix element between two Wilson states is 

j^,^' = j\li:,^^Mm"\lY.,^cj^,m" (5) 

This means that treatment of the Friedel-Anderson and the Kondo im- 
purity can be directly combined, for example, with density functional theory 
calculations of the s-electron band, the Coulomb and s-d-exchange interac- 
tion. One has just to extract from the band-structure calculation the energies 
and matrix-elements of a relatively small number of states (where A^ is of 
the order 40) 
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4 Application of the Modified Wilson States 
in the Kondo Impurity Problem 

As an example I calculate the Kondo ground state and its singlet-triplet 
excitation energy for three examples and compare the results. For the calcu- 
lation I choose (i) the Wilson spectrum with p (£:)=l/2=const, (ii) a three- 
dimensional free electron band with p [e) cc ^Je and (iii) a tight-binding 
band for an fee lattice. The density of states of the latter is sketched in 
Fig. 3. These calculations are an extension of a recent paper on the Kondo 
effect [5] . In the appendix the authors ansatz for the Kondo solution is briefly 
sketched. It is discussed in detail in [3]. 



0.5- 




Fig. 2: Energy band for an fee lattice in the 
tight binding model. The dark shaded region 
represents the half-filled band. 

4.1 Comparison between the different band-structures. 
4.1.1 Rectangular Wilson band 

Wilson's original rectangular band, extending from (—1 : +1) with po = 
1/2, is used. The singlet-triplet excitation energy is plotted in Fig.3a as 
a function of 1/ {2p^J\ It follows essentially an exponential law: AE'^t oc 
exp [-1/ (2po J)] . 



8 



4.1.2 Free electron band. 



Here the s-electrons are treated as free electrons. However, the s-band is 
defined as full (i.e. cut off) at the energy where the electron density of one 

electron per spin and atom is reached. Then the s-band is occupied with 1/2 
electron for spin up and down. This defines the Fermi energy which can be 
expressed in terms of the atomic volume 



When the energy is normalized by the Fermi energy and the Fermi level 
adjusted to zero the upper band edge hes at 2^/^ — 1 w 0.5874. The density 
of states at the Fermi energy per atom and spin is po = 3/4. The singlet- 
triplet excitation energy is plotted in Fig.2 as stars. The (almost) straight 
line lies below the results of the rectangular Wilson band but is parallel to 
it. 

4.1.3 Three dimensional tight binding 

For a face centered cubic metal the tight binding energy is given by 



where Eq is a constant energy, t is the coupling between nearest neigh- 
bors, i.e., the hopping matrix element, a is the cubic lattice constant and 

wave number of the electron. By calculating the en- 
ergy for all k-states in the first Brillouin zone (for a crystal of roughly 10^ 
atoms or k-points in the IBZ) one obtains a density of state which, after 
some smoothing, is shown in Fig.2. The actual values of Eq and t cancel out 
through the normalization with the Fermi energy. The density of states per 
spin and atom at the Fermi energy is po = 1-36 and is much larger than in 
the rectangular Wilson band. The singlet-triplet excitation energy is plotted 
in Fig.3 as triangles. The (almost) straight line lies below the results of the 
rectangular Wilson and the free electron band but is parallel to both. 

All three curves follow essentially the exponential law !S.EstjD cc exp [—1/ (2Jpo)]- 
For the half-band width D the value D = 1 is used in all three curves. 




e (k) ^Eo-At 



cos [-^kr^aj cos (^^kyOj 
+ cos {^kyO) cos {^kzoj + cos {■^kzoj cos [-^kxO) 
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Fig. 3a: The logarithm of the singlet-triplet excitation energy 
In [AEst/D) is plotted as a function of 1/ (2Jpo) for three 
different spectra with different density of states po at the 
Fermi energy. For the half-band width D the value of one 
is used. 

The fact that the three curves in Fig. 3a are essentially parallel to each 
other demonstrates that the density of states po at the Fermi energy is the 
appropriate parameter because otherwise the different curves would form a 
finite angle. On the other hand the fact that the curves for the free electrons 
and the tight binding are shifted with respect to the rectangular band can be 
described by an effective half-band width D* which is different from one. The 
shift is just equal to In (-D*). After adjusting an effective half-band width D* 
the plots of In {AEst/D*) yield a universal curve for the three cases as shown 
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in Fig.3b. 




Fig. 3b: The same plot as in Fig.Sa but using an effective 
lialf-band widtli D*. 

In table I the effective values for the half-band width are collected in 
column 3. In column 4 the effective band width 2D* is multiplied with Pq. 
The result is rather surprising. It yields the value one in all three cases. 

The Kendo energy is essentially given by 

The density of states po appears in denominator of the exponent while the 
half-band width D* is just a pre-factor to the exponential function. From the 
numerical calculation follows that for the very different bands the product 
of Po and 2D* is constant and equal to the number of states (per spin and 
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atom) in the whole band which is one. 



Spectrum 


Po 


D* 


2D>o 


Wilson spectrum 


0.5 


1 


1 


Free Electrons 


0.75 


0.66 


0.98 


Tight Binding 


1.36 


.37 


1.01 



Table I: The density of states po at the Fermi level, 
the effective half-band width D* obtained from the 
Kondo calculation and the product 2D*pQ for bands 
with very different band structure. 

5 Conclusion 

The Wilson states are generalized to electron bands with energy dependent 
density of states and s-d-interaction. These state represent a finite energy 
range and possess the full interaction with the impurity. The Wilson states 
permit us to introduce band-structure data into numerical calculations of the 
Priedel- Anderson impurity and the Kondo impurity. This is demonstrated in 
three examples. These calculations show that the properties of the singlet- 
and triplet state of the Kondo impurity for an arbitrary band can be deduced 
from the results for a rectangular Wilson band and the actual density of states 
at the Fermi level. 

Acknowledgement: The research was supported by the National Science 
Foundation DMR-0439810. 
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A Appendix 

A.l Friedel- Anderson impurity 

With the Wilson bases the Priedel- Anderson Hamiltonian has the form 

N N 

Hfa = ^{'^Sucl^c^a+Eddlda+'^Vsd{iy)[dlc^a+cl^da]}+Und+nd-- (6) 

The author introduced a magnetic solution ^ms of the form 

*M5 = [Aa*o_^al+^ + Bdlal^^ + Cal_^d*^ + Ddld*^] \0a+^0a-i) 

The states aj^ and aj_ arc optimized localized states formed from the Wilson 
states. They determine the full bases |a|_,_| and |a|_|. The many-electron 

state |0a+-|-0a-j.) is given by YYi^i a*_,_^ 11^=1^ ^i-i^o where $o is the vacuum 
state. 

This solution has a considerably lower energy than the mean-field solu- 
tion and yields lower magnetic moments than the mean-field solution. The 
threshold of the Coulomb interaction U to form a magnetic moment is almost 
twice as large as in mean-field theory [?]. 

However, the real ground state is a symmetric state. It is obtained by 
reversing all spins in and superimposing the two solutions. Then one 
obtains 

= *M5 (Ti) T *M5 (it) (7) 

= [Aa*_^a*^^ + Bd*_yal^^ + Cal_yd\^ + Dd*_^d\^] \Qa+}^a-i) 
=F \Aal_^al^^ + Bd*_^a*Q^^ + Ca*Q_^d\^ + Dd*_^d\^ |0„_|0„+;) 

The final solution is obtained by minimizing the ground-state energy by 
optimizing the two localized states and al_. The construction of the 

localized states a^j^ and the full basis |c^l±| ^ ^ optimization are 
described in appendix A3. 

This ground state is determined by the components of the two localized 
states a\_^_ and a\_. Since we generally use 40 Wilson states this requires 80 
components. The resulting ground-state energy and the occupation of the d- 
states is of the same accuracy as a large (1/iV) expansion up to second order 
which requires the calculation of more than 10^ self-consistent amplitudes. 
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A. 2 Kondo wave-function 



The Kondo impurity has a magnetic moment. Therefore either the spin- up 
or spin-down d-state is occupied. Zero or double occupancy of the d-states 
is forbidden. Then the Coulomb interaction transforms into an exchange 
interaction of the form 



Hex 'VaJ 



(8) 



5+M/j(0)vi/^ (0) + ^_vI;T (0) vl/^ (0) 
+5,(vI/|(0)vI/^(0)-vI/|(0)vl/i (0) 

where 5+, Sz are the spin operators of the impurity with spin S = 1/2 
and \E'| (0) and ^'j (0) represent field operators. 

In ref. |2] the approximate solution for the Kondo impurity was intro- 
duced. It is given by equ. . 

(9) 

It uses two sets a bases, 1*^^"} ^^^'"^ orthogonal transformed 

versions of a Wilson basis {cj,}. Again the many-electron state |0a+|0a_|) is 
given by 11^=1^ YYi=i ^i~i^o- The first part of the ansatz consists of the 
product of two n-electron state lOa+^Oa-i) multiplied with the two-particle 
state (-BaQ^|(i* -|- Cd*aQ_^^ . If one reverses all spins in this term one obtains 
the second term (after rearranging the operators). This solution is derived 
from the solution of the Friedel-Anderson impurity in the limit where zero 
and double occupancy of the d-level is zero (for infinite Coulomb potential 
U and Ed = -U/2). 

The final solution is obtained by minimizing the ground-state energy by 
optimizing the two localized states aQ_^_ and al_. The construction of the 

localized states aj_|_ and the full basis |al±| as well as the optimization are 

described in appendix A3. 

If one sets C = C and B = B one obtains the singlet ground state. 
If one permits a free variation of all four coefficients one obtains besides 
the singlet ground state also the excited triplet state with C = —C and 
B = —B. We denote this state as the unrelaxed triplet state. One obtains 
the relaxed triplet state when the opposite signs of the coefficients are locked 
{C = —C and B = —B). Then the minimization of the energy yields the 
triplet state with the lowest energy. In present paper the unrelaxed singlet- 
triplet excitation is plotted for the different bands. 
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A. 3 Construction of the Basis 

For the construction of the state and the rest of basis a| one starts with 
the s-band electrons {c|,} which consist of N states (for example Wilson's 
states). The d'^-state is ignored for the moment. 



• In step (1) one forms a normalized state % out of the s-states with: 



N 



4 = E«°4 



(10) 



u=l 



The coefficients can be arbitrary at first. One reasonable choice is q;° = 

• In step (2) (A'" — 1) new basis states a| (l<i<A'" — 1) are formed 
which are normalized and orthogonal to each other and to Cq. 

• In step (3) the s-band Hamiltonian Hq is constructed in this new basis. 
One puts the state at the top so that its matrix elements are i^oi 
and HiQ. 

• In step (4) the {N — l)-sub Hamiltonian which docs not contain the 
state is diagonalized. The resulting Hamilton matrix for the s-band 
then has the form 



m 

Vfr{l) 



Vfril) 

m 





\VfriN-l) 



Vfr{2) 



m 





Vfr{N 






1)\ 



(11) 



E{N-1) J 



The creation operators of the new basis are given by a new set of 



(0 < i < iV — 1). Again the aj can be expressed in terms of the s-states; 
a] = J2u=i ^u'^i- After the state al is constructed the other states aj 
are uniquely determined. The additional s-d hopping Hamiltonian can 
be expressed in the terms of the new basis, and one obtains the Priedel 
Hamiltonian as given in eq. (??)• The state "^ss is formed, and the 
energy expectation value (of the full Hamiltonian) is calculated. 
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• In the final step (5) the state % is rotated in the A^- dimensional Hilbert 



space. In each cycle the state is rotated in the ( al^a] ) plane by an 



angle Oi^ for 1 < io < N — 1. Each rotation by 6ig yields a new oq 



The rotation leaves the whole basis jag, a] | orthonormal. Step (4), the 

diagonalization of the {N — l)-sub Hamiltonian, is now much quicker 
because the (A^ — l)-sub-Hamiltonian is already diagonal with the ex- 
ception of the io- row and the io-column . For each rotation plane 

(^al,al^^ the optimal with the lowest energy expectation value is 
determined. This cycle is repeated until one reaches the absolute min- 
imum of the energy expectation value. In the example of the Friedel 
resonance Hamiltonian this energy agrees numerically with an accuracy 
of 10~^^ with the exact ground-state energy of the Friedel Hamiltonian 
|3]. For the Friedel- Anderson impurity the procedure is stopped when 
the expectation value changes by less than 10~^° during a full cycle. 




ao — Oq cos 6ig + a] sin 6., 
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